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Abstract 

Atmospheric neutrinos are produced in air showers, when cosmic ray primaries hit the Earth's atmosphere and 
interact hadronically. The conventional neutrino flux, which dominates the neutrino data measured in the GeV to TeV 
range by neutrino telescopes, is produced by the decay of charged pions and kaons. Prompt atmospheric neutrinos are 
produced by the decay of heavier mesons typically containing a charm quark [2EH3]. Their production is strongly 
suppressed, but they are expected to exhibit a harder energy spectrum. Hence, they could dominate the atmospheric 
neutrino flux at energies above ~ 100 TeV. Such a prompt atmospheric flux component has not yet been observed. 
Therefore, it is an interesting signal in a diffuse neutrino search, but also a background in the search for a diffuse 
astrophysical neutrino flux. The sensitivity of diffuse neutrino searches with the IceCube Neutrino observatory has 
reached the level of theoretical expectations of prompt neutrino fluxes, and recent results are presented. 
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Searches for diffuse neutrino fluxes have the goal of 
identifying astrophysical and prompt atmospheric neu- 
trinos in the measured data. The signatures of a dif- 
fuse prompt and extragalactic neutrino signal are their 
energies. The primary background of conventional at- 
mospheric neutrinos shows an energy spectrum, which 
is about one power in energy steeper than the spec- 
trum of primary cosmic rays. The reason is the energy 
dependent competition between the decay and interac- 
tion of pions and kaons. At high energies the decay 
length increases, which makes it more likely that the 
mesons interact instead of producing muons and neu- 
trinos. Prompt mesons have short lifetimes and decay 
into neutrinos independent of their energy and arrival 
direction. Thus, their energy spectrum follows the spec- 
trum of primary cosmic rays. The energy spectrum of 
astrophysical neutrinos is expected to follow the pro- 
duction spectrum at the cosmic ray accelerator, which is 
even harder than the cosmic ray primary spectrum ob- 
served at Earth. Figure [T] illustrates the expected spec- 
trum for prompt and astrophysical neutrinos compared 



to the conventional atmospheric neutrino background. 
A diffuse neutrino analysis aims for a precise measure- 
ment of the observed neutrino energy distribution in or- 
der to search for the characteristic transition to a harder 
spectrum at very high energies. 

The calculation of the prompt atmospheric flux re- 
quires the knowledge of the differential cross sections 
for gg — > cc and ^-fusion, and therefore the knowl- 
edge of parton distribution functions at small values of 
Bjorken x. These are very difficult to measure with col- 
lider experiments Q. A measurement of prompt at- 
mospheric muons and neutrinos could therefore com- 
plement the knowledge of differential cross sections for 
parton interactions and structure functions in accelera- 
tor physics. This makes the search for a prompt atmo- 
spheric neutrino flux an interesting topic. 

The IceCube Neutrino Observatory is a cubic- 
kilometer neutrino detector installed in the ice at South 
Pole El. It is based on the optical detection of neutri- 
nos, which interact with the ice and produce secondary 
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Figure 1: Distribution of primary neutrino energies for simulated 
V/j + Vn for conventional atmospheric, prompt atmospheric and astro- 
physical neutrinos folded with the detection efficiency of this analysis. 

particles. Such charged particles emit Cherenkov light, 
which is detected by optical sensors. These sensors are 
attached to 86 strings, which hold 60 sensors each, with 
a standard vertical spacing of 17 m between sensors and 
a horizontal distance of 125 m between the strings. The 
data analyzed here were taken between May 2009 and 
May 2010 with a livetime of 348 days, when the Ice- 
Cube detector was still under construction and in its 
59-string configuration. This analysis is based on a se- 
lection of muon neutrino events by their characteristic 
track-like signature. The temporal and spatial hit pat- 
tern of light produced along the muon trajectory and 
recorded by the optical sensors allows a directional re- 
construction of the inital neutrino, and the brightness 
of the event is the basis of an energy reconstruction. 
The main background consists of atmospheric muons, 
which outnumber the conventional atmospheric neutri- 
nos by more than five orders of magnitude at trigger 
level. The muon background is removed, by selecting 
high-quality upward-going tracks. Such tracks can only 
originate from neutrinos, which are not absorbed inside 
the Earth. The final sample consists of 21943 events 
with a very small muon background contamination esti- 
mated from simulation to be less than 0.2%. 

Figure [2] shows the reconstructed muon energy loss 
distribution measured with IceCube, which is correlated 
with the initial neutrino energy. The data are analyzed 
in a two-dimensional likelihood analysis, evaluating the 
full information from reconstructed neutrino-induced 
muon energies and arrival directions. The expected 
zenith angle distribution of prompt atmospheric and as- 
trophysical neutrinos is isotropic, while the expected 
distribution of arrival directions of conventional atmo- 




Figure 2: Reconstructed energy loss of + events measured with 
the IceCube detector in its 59-string configuration. The red solid line 
shows the best fit result for an extragalactic diffuse neutrino flux, the 
red dashed line shows the upper limit obtained from this analysis at 
90% CL. The orange dashed line shows the upper limit on a prompt 
atmospheric neutrino component. 

spheric neutrinos shows characteristic features. This be- 
havior is modified by an energy dependent angular de- 
tector acceptance, which increases towards the horizon 
for high-energy events due to absorption in the Earth. 
The use of this angular information increases the sensi- 
tivity of the analysis. The fit is able to identify contribu- 
tions from prompt atmospheric and astrophysical neu- 
trinos in the data sample by comparing it with the ex- 
pected distributions for all three neutrino components. 
The challenge of this analysis is systematic uncertain- 
ties, which distort the expected energy distributions and 
could hide or fake a signal. Such uncertainties are iden- 
tified in the modeling of the detector acceptance as well 
as the theoretical prediction of atmospheric neutrino 
spectra. These uncertainties have been parameterized 
and implemented as nuisance parameters into the fitting 
method, in order to quantify their effects on the analyzed 
distributions in an unbiased way. A particular improve- 
ment in understanding the background expectation of 
atmospheric neutrinos was the inclusion of the cosmic 
ray knee into the neutrino flux calculations 0. 

The measured energy spectrum in Fig. [2] largely 
agrees with the expectation of conventional atmospheric 
neutrinos, which is shown for the standard prediction 
from Honda et al 1 6 1 and for the spectrum modified by 
the best fit nuisance parameters. The high-energy tail 
of the distribution shows a different slope appearing as 
a small excess of events, which is at the level of 1 .8 a 
compared to the expectation of a pure conventional at- 
mospheric neutrino sample. 
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Figure 3: Predictions for v„ + prompt atmospheric neutrino fluxes 
(T] |2j [3j in comparison to the expected flux of conventional atmo- 
spheric neutrinos 1 6 1 . The red shaped area marks the theoretical un- 
certainty on the prediction by Enberg et al. This analysis produces 
individual upper limits valid up to 250 TeV, which is the end of the 
sensitive energy range defined by a worsening of the analysis sensi- 
tivity by 5%. The limits on each of the models are similar in nor- 
malization but follow the slightly different shapes of the models. All 
limits lie within the grey shaded area, which marks the region between 
the lowest and highest flux limit. 



The best fit for prompt atmospheric neutrinos is 
zero, while the best fit for an extragalactic neutrino 
flux is non-zero but non-significant. The upper limit 
at 90% confidence level derived from this analysis on 
the prompt atmospheric neutrino flux is a factor of 4.6 
larger than the expected prompt flux calculated by En- 
berg et al H], which has been modified for an improved 
parameterization of the primary cosmic ray spectrum 
and composition Q. Figure [3] shows the limits on 
several prompt neutrino flux predictions in comparison 
to prompt flux expectations. These limits are slightly 
below the prediction by Bugaev et al. 0, but other 
prompt atmospheric neutrino flux predictions are not yet 
in reach with the current sensitivity. 

The preliminary upper limit derived on a generic as- 
trophysical E~ 2 power law spectrum is E 1 d<S>jdE = 
1.4 • 1CT 8 GeVcrrT 2 s _I sr _1 , which is slightly above the 
Waxman-Bahcall bound. The limit on a diffuse astro- 
physical flux is presented in Fig. [4] 

The completed IceCube detector will provide much 
higher statistics than this data sample and an expan- 
sion to the higher-energy regime. Events with cascade 
signatures allow a very precise energy reconstruction 
and smaller background rates from conventional atmo- 
spheric muons and neutrinos. In a combined analysis, 
IceCube will be able to challenge prompt neutrino flux 
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Figure 4: Limit on an extragalactic diffuse neutrino flux from this 
analysis in comparison to limits from previous experiments |7][8). The 
red solid line marks the Waxman-Bahcall upper limit as a benchmark 
for extragalactic diffuse neutrino expectations 1 9] 1 101 - The dashed 
lines show flux predictions for extragalactic neutrino fluxes for differ- 
ent source classes II 1 II 1 2llT3l [T4lfT5l ■ 

predictions as well as astrophysical flux models within 
several years of operation. 

References 

[1] R. Enberg, M. H. Reno, I. Sarcevic, Phys. Rev. D 78 (2008) 
043005. 

[2] A. D. Martin, M. G. Ryskin, A. M. Stasto, Acta Phys.Polon.B 

34 (2003) 3273-3304. 
[3] E. Bugaev, V. Naumov, S. Sinegovsky, E. Zaslavskaya, II Nuovo 

Cimento 12C 1 (1989)41. 
[4] A. Achterberg et al. (IceCube Collaboration), Astroparticle 

Physics 26 (2006) 155-173. 
[5] T. K. Gaisser, arXiv: 1111 .6675. 

[6] M. Honda, T. Kajita, K. Kasahara, S. Midorikawa, T. Sanuki, 

Phys.Rev.D 75 (2007) 043006. 
[7] J. A. Anguilar et al. (ANTARES Collaboration), Phys.Lett.B 

696 (2011) 16-22. 
[8] A. Achterberg et al. (IceCube Collaboration), Physical Review 

D 76 (2007) 042008. 
[9] E. Waxman, J. Bahcall, Phys. Rev. D 59 (1998) 023002. 
[10] E. Waxman, To be published in Astronomy at the Frontiers of 

Science. 

[11] K. Mannheim, Astropart. Phys. 3 (1995) 295. 

[12] J. K. Becker, P. L. Biermann, W. Rhode, Astropart.Phys. 23 
(2005) 355-368. 

[13] F. Stecker, Phys.Rev.D 72 (2005) 107301. 

[14] A. Muecke, R. Protheroe, R. Engel, J. Rachen, T. Stanev, As- 
tropart.Phys. 18 (2003) 593-613. 

[15] S. Razzaque, P. Meszaros, E. Waxman, Phys.Rev. D 68 (2003) 
083001. 



